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ABSTRACT 

In the present research , microwave induced hyperthermia technique was used as a solution , for the treatment 
of breast tumor of volume 561300 mm 3 . Techniques like chemotherapy , radiotherapy or surgery are generally used for 
the treatment of tumor. The present work describes the simulation carried on actual tumor images obtained from MRI 
scans of the breast. An optimized heat input (1W - 4W) was selected for constant frequency with respect to transient 
response Analysis is done on 2D finite element generated pathological human tissue containing tumour. The operating 
frequency of the antenna is 2.4GHz. Heat is generated using microwave induced hyperthermia and simple methodology 
of finding effective model parameters for the optimal interstitial microwave hyperthermia treatment is determined. 
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INTRODUCTION 

The term hyperthermia refers to either a phenomenal high fever or the treatment of an infection, by the 
infusion of an abroad protein utilizing heat. In Hyperthermia temperature of a piece or possibly the entire body 
raises above the typical degree for a characterized timeframe. The level of temperature, height related to 
hyperthermia is around a couple of degrees above the ordinary temperature (41-45° C) (Habash et al. 2006). 
Hyperthermia used with chemotherapy or radiation therapy elevates the body temperature above the ordinary 
temperature (41-45°C) to cause damage to cancer cells [1]. Different types of energy may be used to apply heat, 
in hyperthermia including microwave, laser, magnetic field and ultrasound [2]. In the Hyperthermia treatment 
heating of the tumour is always done under controlled conditions. Since, healthy tissue can absorb energy, 
causing damage with the application of microwave, laser and ultrasound. 

Treatment by Hyperthermia 

The best remedial use of heating, that have physiological most essential hyperthermia is the treatment of 
hurts, agonies, stresses, and sprains by means of use of temperature beneath 41°C for 60 minutes or thereabouts 
and utilize physiological parts of expanding blood flow and metabolic rates (Roemer 1999). 

For tumour treatment purposes, there are reports that malignant cellular material is usually more 
sensitive to temperature than are their normal counterparts 

According to a study by American cancer institute 55% cases of cancer are treated with Surgery, 18% 
with chemotherapy, 11% with radiotherapy and the remaining 16% using other techniques. 
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Techniques used for Hyperthermia Treatment of Cancer 
Magnetic Heating 

The National Cancer Institute classifies hyperthermia into three different methods on which research is being 
done. They are 

• Local hyperthermia 

• Regional hyperthermia 

• Whole body hyperthermia. 

Local hyperthermia is a technique in which a local area such as tumor is treated by using various methods that deliver 
heat to the infected part. Heat may be applied using different forms of energy which can be radio frequency, ultrasound, 
and microwave. To treat the tumors just below the skin external approaches are used and to treat the tumors within or near 
the body cavities Intra luminal or endocavitary methods are used, whereas, tumors deep within the body are treated by 
interstitial methods. 

In regional hyperthermia large area such as body cavity, the organs are heated using deep tissue, regional perfusion 
and continuous hyper thermic peritoneal perfusion. 

And Whole body hyperthermia is treatment of cancer which has spread throughout the body. Since the body can 
tolerate the temperature rise up to 45°C, therefore till this temperature is attained normal tissues are not damaged. But some 
tissue may be affected due to difference in temperature rise at various spots which may cause burns, blisters, pain, 
discomfort etc. [3-4] 

Microwave Hyperthermia 

It is one of the promising techniques which are being used on hundreds of patients, those who are suffering from 
either breast cancer or prostate cancer. Microwave energy is effective in heating cancerous tumors because tumors have 
high water content. One or more microwave antenna can be used to treat tumor that depend on tumor size and location in 
the body. Microwave induced hyperthermia treatment utilizes a single waveguide microwave antenna working at 434,915 
and2450MHZ [5-6]. 

Background of Work Done and Methodology Adopted 

The breast cancer ablation by radiofrequency has been associated with the capacity of concentrating the 
electromagnetic energy in the tumour area [7]. Radiofrequency ablation in the recent years reached importance for the 
treatment of cancerous tumor. Breast cancer is believed to have resulted from the interaction of genetic factors, with other 
environmental factors. The available techniques for breast cancer treatment have strong side effects. The most common 
used treatment is a mastectomy, an approach which comes down to removing the entire breast and the other is a 
lumpectomy, in which only the one part of the breast is removed. Other techniques used are such as radiation therapy and 
chemotherapy that cause the cellular damage all these are ablative techniques. In these treatments the tissue temperature is 
raised in order to kill the malignant tissue. These types of treatments have a lot of potential, since it can offer a non- 
invasive treatment with low side effects [8]. In hyperthermia temperature is raised above 42-43°C for the sufficient period 
of time. The challenge of achieving the temperature with non-invasive electromagnetic hyperthermia treatment is to focus 
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electromagnetic (EM) power in the cancerous tissue while avoiding the normal tissue. In microwave frequency, energy 
couples into the tissue through antennas that emits the microwaves [9]. The challenge while working with the microwave 
ablation is to avoid heating the healthy tissue while heating the tumor. While the relatively high conductivity of tumor 
tissue increases the local heating potential and thus high temperature difference can be obtained. The hot spots in the 
healthy tissue areas results in undesired side effects such as burns, blisters and extra pain. More flexible approach, focusing 
the EM energy is the use of the phase array system. [10] .A 2D finite element analysis for the treatment of pathological 
human tissues containing tumors is presented. The analysis compares the model of coaxial antenna with one, two and three 
air slots. The model is based on the TM mode coupled with Penne’s equation under a transient state condition. The 
operating frequency of the antenna using is 2.45GHZ [11]. A zero order mode resonator (ZOR) meta material (MTM) 
structure for microwave thermotherapy to create and transmit a plane EM wave into treated natural tissue, that can deliver 
a homogeneous SAR circulation in the plane parallel to tool gap, is generated. ZOR principle is simulated in COMSOL 
software that shows a very good SAR homogeneity and achieves a penetration depth close to that of EM plane wave [12]. 
A method of microwave hyperthermia by utilizing thin microwave antenna located in human tissue is described. For the 
numerical simulation finite element method (FEM) is used. 

STRUCTURE OF ANTENNA 

The antenna consists of a thin coaxial cable with a slot measuring 1 mm cut on the outer conductor 5mm for the 
tip. The antenna is enclosed in a catheter made of poly tetra fluoro ethylene) and operates at the frequency of 2.4GHz, 
which is used in microwave coagulation therapy [13]. 


Table 1: Dimension of Coaxial Antenna 


PROPERTIES 

VALUE 

Diameter of catheter 

1.79mm 

Diameter of the central conductor 

0.29mm 

Inner diameter of the outer conductor 

0.94mm 

Outer diameter of the outer conductor 

1.79mm 


Table 2: Properties of Material 


PROPERTIES 

CONDUCTOR 

(Dielectric) 

CATHETER 

Relative permittivity 

2.03 

2.60 


REAL IMAGE GEOMETRY MODEL 

In FEM analysis, a correction of the model plays an important role in the simulation. In this paper a set of real 
image from the MRI sliced was introduced to enhance the breast model design. The actual images help in working on 
actual volume and size of tumor for better result. The motive is to find out the time in which required temperature can be 
obtained, which is possible only with the actual volume and size of the tumor. MRI images in axial view, sagittal view and 
coronal view were imported and rebuild to obtain breast model, with .stl format. This was done with the help of mimics 
software and through manual segmentation. Later the model is imported into COMSOL multi physics for the simulation. 

MODEL SIMULATION 

As described above, the breast model is developed from a set of MRI images and, then the microwave antenna 
was modelled and inserted into the tumor tissue. The destructive area of the tissue was finally expected to be identified. 
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The microwave antenna is kept in the frequency of 2.4GHz. The heat generated from the electric field in the waves is 
analysed by using COMSOL Multi physics. Modules selected for the simulation are RF module and the bio heat. FEM 
meshing refinement is done by automatic initial mesh of tetrahedral mesh shape. The physics and all the boundary 
conditions are assigned to the domains and are shown above in the table 1. Figure 1: Domain represents the breast tissue, 
whereas the circled domain is the tumor inside the breast model 



Figure 1: MRI Image of the Cancer Patient 



Figure 2: 3D CAD Model of the Cancerous Breast 
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Figure 3: Meshing Refinement in the both the Domains (Breast and the Tumor Region) 



Figure 4: Meshing Refinement of the Cancerous Tumor 


Table 3: Physical Properties of Tissue 


PROPERTIES 

VALUE 

Heat Capacity 

2348 [J/kg.K] 

Density 

911 [kg/m 3 ] 

Thermal Conductivity 

0.21 [W/m.K] 

Frequency factor 

4.436*10 16 [1/s] 

Activation energy 

1.3*10 5 [J/Mol] 


Table 4: Physical Properties of Tumor 


PROPERTIES 

VALUE 

Heat Capacity 

3700 [J/kg.K] 

Density 

1050 [[kg/m 3 ] 

Thermal Conductivity 

0.5 [W/m.K] 

Frequency factor 

9.4*10 1U4 [1/s] 

Activation energy 

660000 [J/Mol] 


DOMAIN AND BOUNDARY EQUATIONS 
ELECTROMAGNETICS 

Coaxial cable caring electromagnetic wave is characterized by transverse electromagnetic fields (TEM). 
Assuming time harmonic fields with complex amplitudes containing the phase information, the appropriate equation is: 
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E = e r £ e i< “ t_kz) 

r 

H = e (p ^e i<m, ' kz) 


p av = Sinner R < E * H)2mdr = 


C A 2 t / router x 
Z rinner' 


Where, 


z 

Direction of propagation 

r,cp and z 

Cylindrical coordinates centered on 
the axis of the coaxial cable 

Pav 

Time averaged power flow in the 
cable. 

Z 

Wave impedance in the dielectric of 
the cable 

r inner and r outer 

Dielectric inner and outer radius 

Q 

Angular frequency 


The propagation constant k, relates to the wavelength in the medium, X as: 

7 2n 

k = T 

In the tissue, the electric field also has a finite axial component, whereas the magnetic field is purely in the 
azimuthal direction. Thus, model the antenna using an axisymmetric transverse magnetic (TM) formulation. 
The wave equation then becomes scalar in Hep: 

V * (Ye r - * H<p) - iV^Hcp = 0 


The boundary conditions for the metallic surfaces are: 


N*E = 0 


The point is modelled using a port boundary condition with a power level set to 10W. This is essentially a first 
order low, reflecting boundary condition with an input field H(p 0 : 

n * VsE - VjUH(p= ^V^Hcpo 

r _ PavZ _ 

Where, H»„ = wr/rmner j 

r 

Applying this boundary condition to all boundries 

HEAT TRANSFER 

The bio heat equation describes the time dependent heat transfer problem as: 

P Cpf + V (-kVJ) = p b C b co b (T b -T) + Q met + Qext 
Where, 
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k 

Tissue thermal conductivity (W/m.k) 

Pb 

Blood density (kg/m 3 ) 

c b 

Blood specific heat capacity 

C0 b 

Blood perfusion rate (1/s) 

T b 

Arterial blood temperature (K) 

Qmet 

Heat source from metabolism (W/m 3 ) 

Qext 

External heat source (W/m 3 ) 


The initial temperature T b is applied to all domains. The heat source from the metabolism is neglected in this 
model. The external heat source is equal to the resistive heat generated by the electromagnetic field: 

Qext = '/2 Re[(o - jcas)E.E] 


Table 5: Parameters of Blood 


PARAMETERS 

VALUE 

Density of blood 

1000[kg/nr’] 

Heat capacity 

4180[J/kg.K] 

Blood perfusion rate 

0.0064[l/s] 

Arterial blood temperature 

37° 

Initial and boundary temperature 

37° 


In addition to the heat transfer, the tissue damage integral is also computed. This gives an idea about the degree of 
tissue injury a during the process, based on the Arrhenius equation: 

da . , dE x 

— = A exp(-—) 

dt r RT' 

Where, A is the frequency factor (1/s) and De is the activation energy for irreversible damage reaction (J/mol). 

SPECIFIC ABSORPTION RATE DISTRIBUTION (SAR) 


SAR can be defined as the energy absorbed by living tissue when the body is exposed to electromagnetic field at 
radio-frequency range can be measured by the specific absorption rate. This indicates the absorption power per mass of 
living tissue or watts per kilogram (W/kg). The heating ability of microwave antenna is studied with the help of SAR 
distribution. Heat generated by an electric field in living tissue is given by: 

c t*E a 2 

SAR = —— 

p 


Where, 


P 

-55 - 

Density of tissue (kg/m ) 

a 

Conductivity of tissue (S/m) 

E 

Electric field (V/m) 


Heat generated by the antenna as SAR can be applied as an external heat source in bio heat. Therefore, bio heat 
equation can be written as: 

pCp ^ + V (-AVT) = h(T b - T) + p.SAR 


RESULTS 

Figure 5 and Figure 6, shows the temperature distribution for the duration of 10 min. and 20 min., respectively in 
the tumor, at the frequency range of 2.4GHz. The peak value of the temperature is around the slot and its highest value 
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near the microwave antenna slot. Maximum temperature of 52°C is attained which is optimum for the treatment of the 
tumor of microwave induced hyperthermia. Temperature in a tumor model rises to 52°C within 10 minutes at power of 3W 
and it is observed after a certain period of time temperature remains constant at the same point. Therefore, the temperature 
after 20 minutes remains 52°C. A temperature higher than 52°C can be harmful to the body where it may affect the healthy 
tissues. Temperature distribution on the surface of tumor at three different points during the duration of 10 minutes is 
shown in Figure 7 and these points are located at random location according to the Figure 17from the tip of the probe. 
Similarly, Temperature distribution on the surface of tumor at three different points during the duration of 20 minutes is 
shown in Figure 8. Table 6shows the time required for attaining particular temperature and the power corresponding to it, 
and represents temperature at a point when the power supply of the antenna is varied from 1W to 5W and the duration of 
time which is used for the heating of the tumor are 10 min, 20min, and 30 min, respectively. It is clear from the readings 
that, the temperature rise in the tumor is seen up to 10 min and after that, the temperature inside the tumor is constant and 
does not vary in 20 & 30, with an increase in time. The same is seen for the other power supply values that’s as we 
increase the power of the antenna the temperature increase up to 10 minutes and after that the temperature becomes 
constant. Figure 9 shows the SAR distribution in the breast model. As it can be seen that the sudden peak in the graph 
shows the sudden increase in the value of SAR, this sudden peak can be seen in the area closer to the tip of the antenna and 
as the value of SAR is evaluated away from the antenna it come out to be less as which can be observed in the plot. The 
power of the antenna is kept, such as the there is no effect of electromagnetic wave on the other parts of the body other 
than the part that is diagnosed that’s the tumor region. The figure 11, 12, 13,14,15,16 shows the current density, total 
power dissipation, electric displacement field and electric field in the breast and tumor model. The result is computed at the 
breast model whose volume is 561300 mm 3 and the surface area of the breast model is 3600 mm 2 , whereas the volume of 
the tumor region is 1734 mm 3 and the surface area of the tumor region is taken as 877.7 mm 2 . The power required and the 
temperature produced inside the tumor also depends upon the area of the tumor region. 

Table 6: Temperature Rise Corresponding to Power and Time 


POWER (W) 

TEMPERATURE (max.) 

10 min 

20 min 

30 min 

1W 

43°C 

43°C 

43°C 

2W 

48°C 

48°C 

48°C 

3W 

54°C 

54°C 

54°C 

4W 

60°C 

60°C 

60°C 

5W 

65°C 

65 °C 

65°C 
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Figure 5: Temperature Distribution in the Breast Model 



Figure 6: Temperature Distribution in the Tumor Model 


Point Graph: Temperature CdegC) 



Figure 7: Temperature Distribution on the Surface of Tumor at Three 
Different Points during the Duration of 10 Minutes 
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Point Graph: Temperature (degC) 



Time (min) 


Figure 8: Temperature Distribution in the Surface of Tumor at Three 


Different Points during the Duration of 20 Minutes 



Figure 9: SAR Distribution in the Tumor at the Slice View 


Line Graph: emw.Qh/rho_blood (W/kg) 



Figure 10: SAR Distribution Plot at the Surface of the Tumor, the Antenna 
Inserted Inside the Tumor is at a Frequency of 2.4GHz 
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xlO 3 



Figure 11: Current Density in the Breast and Tumor Model in the Slice View 



Figure 12: Total Power Dissipation in the Breast and Tumor Model in the Slice View 


Slice: Electric displacement field, x component (C/m*) 



xlQ" s 



0.5 


0 


■ - -0.5 

l 


Figure 13: Electric Displacement Filed in the x Direction in the 
Breast Tumor Model in the Slice View 


To find out the relation between volume of breast and the power required, new term a is introduced. 
Where a = (volume of breast/power required) /1000 

Table 7: Shows the Value of a at Different Power and Time Required to Attain the Same 


Time (m) 

Temperature (°C) 

l 


1 W 

2W 

3 W 

4 W 

0 

37.00 

37.00 

37.00 

37.00 
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Table 7: contd., 


1 

38.63 

40.26 

41.89 

43.52 

2 

40.14 

43.28 

46.40 

49.56 

3 

40.97 

44.89 

48.87 

52.85 

4 

41.41 

45.75 

50.23 

54.69 

5 

41.67 

46.23 

51.00 

55.71 

6 

41.82 

46.54 

51.47 

56.32 

7 

41.92 

46.76 

51.77 

56.71 

8 

41.99 

46.91 

51.97 

56.98 

9 

42.03 

47.02 

52.11 

57.16 

10 

42.07 

47.10 

52.21 

57.29 

a 

11.46 

5.73 

3.82 

2.865 


Point Graph: Temperature (degC) 



Figure 14: Temperature Rise at Power of lw 


Point Graph: Temperature (degC) 


l3 



Figure 15: Temperature Rise at Power of 2 w 
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Figure 17: Point in the Tumor where Temprature is Measured at 3W Power 



Figure 18: Electric Field in the Breast and Tumor Model in the Slice View 


CONCLUSION AND DISCUSSIONS 


The results present a simple methodology of finding best model parameters for the optimal interstitial microwave 
hyperthermia treatment. A rise in temperature of breast model was observed because of the RF heating. 

We are able to see a rise in the temperature in the breast model because of RF heating. The intensity of power of 
the antenna is kept between 1W to 3W and the maximum temperature attained by the breast model is between 43°Cto 54°C 
and respectively the time duration under which the temperature increase was studied was 10, 20, 30 minutes respectively. 
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The temperature attained is sufficient for eliminating the cancerous cells. The study of the SAR (Specific Absorption Rate) 
is also done. 

Foster, et. al., recommended uniform semi-infinite homogeneous and uniform homogeneous sphere plane of the 
tissue. SAR calculation and the resulting temperature rise were projected by the finite element solution of the bioheat 
equation. The effect of the distance between head model and the RF source of SAR distribution is inspected. It was 
determined that the temperature elevation reached its steady state value after several minutes [14] .Abdallah et al. 
Investigated the SAR effect in the frequency range of 100MHz to 300MHz in perpendicular polarized plane waves. The 
result was computed using the analogy between uniform transmission lines where each layer is replaced with the 
equivalent transmission line of the same impedance [15]. 

We are able to study the hyperthermia treatment of cancer by converting the MRI images into the 3d model and 
an attempt has been made to study the relationship between the temperature and its effect on the tumor region. 
The observation is done at the different points which are 1mm, 2mm, and 3mm away from the tip of the antenna. As 
studied in literature for the effect of the temperature of the cancerous tissue. The effect of the temperature was studied for 
the different duration that is for 10, 20, and 30 minutes at the different power supply that is from 1W to the 5W. 

Both electromagnetic and ultrasound cylindrical phased array can be used for the hyperthermia treatment of tumor 
within the breast model. The US designed applicator has a smaller focus size and larger penetration depth than that of the 
existing hyper collar. Therefore, US applicator is more appropriate for the smaller and the deep located tumor regions, 
whereas for the larger size tumor hyper collar is preferred and it performs well. Therefore, the further research is required 
to investigate for the heating profiles. 

From the temperature profile distribution, it can be seen that focusing point of the electromagnetic field near the 
tip of the antenna which is responsible for the hot spot in the temperature distribution in the tissue. The material used for 
the antenna is dielectric and catheter. The frequency at which the power is supplied to the antenna is taken as 2.45GHz. 

Our study will continue with the investigation of non-homogenous heating. For example, a tumor as a sub domain 
embedded in the normal tissue, or the insertion of ferromagnetic particles in the targeted tissue, in order to better 
concentrate the radiation energy into it. Also, the study of multiple pin applicator, such as the array of implanted probes 
investigated for the extension of the conclusions adopted. 
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